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ABSTRACT 

\W 

A sudden increase of the particle intensity and of the magnetic 
field strength on September 25, 1961 at 17:31 U.T. was detected by the 
ion- electron detector and the magnetometer on board Explorer XII out- 
side the earth's magnetosphere. From the physical extent of the region 
and the energy of the trapped protons the temporal nature of the event 
was established. Similarities between particle fluxes as well as be- 
tween magnetic fields during the increase and within the outer belt 
make it probable that the increase represents an impulse-like exten- 
sion of at least a portion of the outer belt. A deduced lower limit on 
the outward radial velocity of the extension is 156 km/sec. Crude cal- 
culations show that either the proton density near the outer edge is 
less than 100 cm -3 or the temperature is greater than 15 ev. 
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and 
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INTRODUCTION 

The Explorer XII satellite was well suited for the study of the outer 
boundary of the trapping region. With an apogee of over 13 R e (earth* s 
radii) it usually penetrated the boundary on the sunward side twice 
during each orbit. The position of the boundary ranged from about 8 R e 
to beyond the apogee of the satellite with an average value of some- 
where about 10 R e (Cahill and Amazeen, 1963). 

Data from the ion- electron detector flown by Davis and Williamson 
and from the magnetometer flown by Cahill on board Explorer XII has 
yielded more than 100 transitions of the outer boundary. For our pur- 
poses a transition shall be defined to occur when on an inbound pass the 
output current of the photomultiplier, measured in the mode which is 
most sensitive to the incident energy flux, rises significantly above the 
dark current (usually at least by an order of magnitude) in a time inter- 
val of the order of minutes corresponding to a radial distance of about 
one hundred kilometers covered by the satellite. The reverse is of 
course true for an outbound pass (Davis and Williamson, 1963). Simi- 
larly, a transition of the boundary is indicated by the magnetometer 
when there is usually an abrupt change in either both the magnitude and 
direction of the magnetic field, or just a change in the direction of the 
magnetic field. Several typical transitions are shown in Figures la 
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Figure la-Three samples of Explorer XII transitions of the outer boundary of the trapping region 
as seen by the ion-electron detector. The 8th dynode current is a measure of the electron inten- 
sity for electrons with 20 kev < E < 1 Mev. 
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Figure Ib-Magnetic field directions and intensity for the same transitions as in Figure la. The 
strength is measured in gammas and the directions are given in the satellite centered spherical 
coordinates a and \p. 


and lb. As can be seen from Figure la, at times the satellite passes 
through regions of enhanced particle intensity which are separated from 
the boundary by stretches where the particle intensity falls below the 
threshold of detectability. Unfortunately, in most cases it is impossible 






















to separate temporal and spatial phenomena on the basis of only particle 
data. One could argue equally well that the intensity represent actually 
a part of the outer belt that for some reason is separated from the rest 
of the belt by particle free regions, or that they represent an expansion 
and contraction of the trapping region in such a way that the satellite 
crosses the boundary several times during one orbit. Figure lb shows 
the magnitude and direction of the magnetic field during the same tran- 
sitions. The directions of the magnetic field are given in satellite 
centered spherical coordinates (Cahill and Amazeen 1963). It should be 
noted that in analogy to particle intensities during the transition of the 
boundary the magnetic field sometimes shows the same direction and 
magnitude in regions separated from the magnetosphere as were ob- 
served in the trapping region. Using magnetic field data alone one 
again cannot distinguish between temporal and spatial phenomena. 

In the analysis of the data thus far, one case was found in which a 
combination of magnetic field data and particle data can resolve this 
question. 

The present paper shall be concerned with a sudden increase of the 
particle flux on September 25, 1961 at 17:31 U.T. at a radial distance of 
77,465 km as seen by the ion-electron detector and changes in the mag- 
netic field as seen by the magnetometer on board Explorer XII. 


INSTRUMENTATION 

A brief review of some of the features of the ion-electron detector 
pertinent to the following discussion will be given here. A detailed 
description of the apparatus will be published elsewhere (Davis and 
Williamson, to be published). 

The detector is a scintillation counter consisting of a photomulti- 
plier tube on the face of which is deposited a 5 mg/cm 2 thick layer of 
crystalline ZnS covered by a 1000 A° thick layer of Al. 

It has a geometric factor of 5.4 x 10 -3 cm 2 ster. A collimator 
allows particles to enter from a viewing cone with a half angle of 1 1 . 

A wheel driven by a stepping motor introduces varying thicknesses 
of Ni absorber between the collimator and the phosphor. Three posi- 
tions on the wheel carry Au discs which scatter the incident flux enter- 
ing through an alternate collimator into the phosphor. 
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There are three modes of operation: 

1. Pulse output from the anode measures protons between the 
energies of 100 kev and 5 Mev. 

2. 8th dynode current measures the total incident energy flux from 
both protons and electrons. 

3. Electrons can be selected preferentially by scattering an inci- 
dent beam of both protons and electrons from an Au disc. In 
this mode of operation the 8th dynode current measures the 
total incident energy flux due to electrons alone. 

In all three modes of operation the lower energy cutoff can be raised 
by introducing varying thicknesses of Ni absorber. 

Throughout the following discussion we shall make the assumption 
that for different absorbers we are measuring the integral proton flux 
above some energy for which the detector efficiency is 50%. The 
situation is more complicated for electrons because of the shape of the 
sensitivity curves and the necessity to use approximation methods to 
obtain electron spectrums. 

The magnetic field data was obtained from the magnetometer flown 
by L. J. Cahill of the University of New Hampshire and its operation 
had been described earlier (Cahill and Amazeen, 1963). 


ANALYSIS 

The inbound pass of September 25, 1961 had a geomagnetic latitude 
of less than 5 for radial distances of more than 40,000 km. Figure 2 
shows intensities of low energy protons and the value of the photo- 
multiplier 8th dynode current obtained with 5/j. inch Ni absorber when 
the detector was looking at particles with local pitch angles of around 
90° . Also shown are the intensity and direction of the magnetic field 
in the satellite coordinate system. 

The electrometer data represents instantaneous values of the 8th 
dynode current sampled at intervals of 83 seconds. To gain some 
physical idea of the electron intensities involved, a crude analysis 
shows that at 19:45 U.T. the differential electron spectrum could be 
approximated by a curve of the form E“ r where E is the electron 
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Figure 2- Inbound pass of September 25, 1961. The upper two curves represent the unidirectional 

proton intensity for particles with 90° local pitch angle and the 8th dynode current corresponding T 

to electrons with E > 20 kev and 90° pitch angle respectively. The lower three curves show the 

strength and directions of the magnetic field. Angles a and i//were previously defined by Cahill 

and Amazeen (1963). The left side is a typical response for the trapping region (Region I). The 

right side (Region II) shows the usual response of the detectors outside the magnetosphere. For • 

protons the threshold of sensitivity is the count rate corresponding to an omnidirectional flux of 

penetrating particles. Region III is the sudden increase in particle and magnetic field intensity 

which is the subject of the present study. 
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energy and y - 3.7. The integral intensity of electrons for E > 20 kev 
is 4 x 10 6 electrons (cm 2 sec ster) -1 . This number is accurate to 
better than a factor of 2. 

The plot of the proton data was made for particles with energies 
E > 140 kev and E > 200 kev. The proton intensities are obtained from 
data accumulated over a time period of 0.3 seconds. During that time 
the satellite turns about its axis through an angle of about 58 . Thus 
the 90° pitch angle intensity is an average value covering also other 
than locally mirroring particles. 

The magnetometer data corresponds to an average of ten succes- 
sive telemetry samples taken at intervals of 50 km. The field strength 
is expressed in gammas, and the direction is indicated in spherical 
satellite centered coordinates mentioned earlier. 

The graph in Figure 2 is divided into three parts: Region I repre- 
sents the proton and current intensities in the outer belt. Region II is 
characteristic of the space outside the trapping region. In this case 
the current output of the detector consists only of the dark current and 
the pulse output merges into the omnidirectional background. Region HI 
shows a sudden increase in the electron and proton intensity at 17:31 U.T. 
which lasted for about one minute. 

Similarly the magnetic field in Region I is the rather constant field 
encountered in the outer belt. Region II is typical of the field outside 
the boundary of the magnetosphere. The field is rather disorganized 
and shows many changes in direction. The dashed line indicates the 
magnetic field in Region III. During the transition of Region HI the mag- 
netic field agreed in magnitude and in direction with the field observed 
in Region I. Before and after the sudden increase the magnetic field had 
an average value of 25 gamma. During the transition of Region III the 
field rose to 61 gamma. Aside from this single increase there is nothing 
unusual happening in the behavior of the field outside the magneto spheric 
boundary. It should be emphasized that in this context the word "region" 
is used only with reference to the graph and does not necessarily imply 
any physical structure in space. The increase is real and not due to 
faulty data processing because it appears on two completely independent 
data tapes acquired at Fort Myers, Florida and at South Point, Hawaii. 
This increase in particle and field intensity will from now on be re- 
ferred to as Region III. 
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While the satellite was in Region III it was possible to obtain 15 
data points for protons with E > 140 kev and 45 data points for protons 
with E > 200 kev. Figures 3 and 4 present the pitch angle distributions 
for these two cutoff energies. In both cases the pitch angles were as- 
signed from the knowledge of the direction of the magnetic field as seen 
by the magnetometer. The vertical bars represent Poisson standard 
deviations only. The uncertainty due to the digitization of the analog 
data tapes is about half as large again. In Figure 3 the horizontal bars 
represent the ranges of pitch angles which were scanned during the 
accumulation period. The pitch angle assigned to each point is the value 
seen at the center of the accumulation period. The points represented 



Figure 3-Pitch angle distribution of protons with E > 140 kev found in Region III. The 
vertical bars represent Poisson standard deviations. The uncertainties due to the digi- 
tization of the analog data are about half as large again. The horizontal bars repre- 
sent the range of angles over which the samples were collected. The assigned p.a. 
corresponds to the center of the accumulation period. 
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Figure 4-Pitch angle distribution of protons with E > 200 kev. The data points were 
lumped together into three groups each containing 15 points. The vertical bars repre- 
sent Poisson deviations. The background count obtained with a thick plug in front of 
the detector aperture is given below. 


by crosses actually have pitch angles which are supplements of the 
indicated pitch angles. Figure 4 is a similar diagram. In this case the 
number of counts per data point was small so that data points were 
lumped together into three groups each containing 15 points. The bars 
again represent Poisson standard deviations. For contrast the omni- 
directional background count obtained with a thick plug in front of the 
detector is plotted below. 

The shape of the pitch angle distributions is quite similar to those 
obtained for trapped particles within the magnetosphere in that the 
maximum intensity is found at 90° local pitch angle (Davis and 
Williamson, 1963). By contrast, solar protons during a sudden com- 
mencement have been observed to have maximum intensities in direc- 
tions other than 90° to the field lines (Hoffman, et al. , 1962). 
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Figure 5 shows the magnetic field intensity and directions during 
the transition of Region III. The vertical bars represent standard de- 
viations of the mean of 16 data points each. The scatter in the original 
data is mostly due to the digitization of the magnetometer analog output 
voltage. The standard deviations increase as the field decreases in 
magnitude, a is the angle between the spin axis of the satellite and the 
magnetic field. \p is the angle between the spin-field plane and the spin- 
sun plane. They are identical to those defined earlier by Cahill and 
Amazeen (1963). Before and after the transition of Region HI the mag- 
nitude of the field was about 25 gamma. During the transition the mag- 
nitude rose to 60 gamma. There was no drastic change in either a or 0. 

Assuming for the moment that Region III is stationary an esti- 
mate of its maximum width can be made by finding the time when the 
particle flux was zero on the step corresponding to the lowest energy 
cutoff just before and after the transition of Region III and multiplying 
it by the velocity of the satellite. This gives the region a maximum 
width of 200 km. A more precise value obtained from the magnetometer 
data is 83 km corresponding to a width of 56 km in the radial direction 
from the center of the earth. It should be pointed out that the angle be- 
tween the direction of the magnetic field and the velocity vector of the 
satellite is about 60 . This means that the angle between the velocity 
vector and the plane of gyration of protons with 90° pitch angle is 30°. 
Therefore, the width of Region III in the plane containing protons with 
90° pitch angles must be close to that of Region III as measured along 
the velocity vector. 

Table I shows the cyclotron radii of protons with several energies 
corresponding to the magnetic field encountered within Region III. It 
is clear from the Table that no 200 kev or 140 kev protons with pitch 
angles close to 90° could be trapped within a stationary Region III where 
the magnetic field strength is about 60 gamma. The efficiency of the 
detector for 310 ev protons is zero. Thus the conclusion must be drawn 
that Region III is a temporal rather than a spatial phenomenon. 

The origin of Region III can only be conjectured. In what follows 
we shall try to produce convincing evidence that Region III represents 
a temporal extension of the outer belt. 

Comparison between Region I and Region HI shows that the intensi- 
ties of both the 8th dynode current and of the protons are comparable 
in several places for particles with 90° pitch angle. The order of 
magnitude fluctuations of the proton flux and the somewhat smaller 
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Table I 


Cyclotron Radii of Protons 


Magnetic field strength 

H = 60 gamma 

Proton energy 

Cyclotron radius 

ZOO kev 


1075 km 

140 kev 


899 km 

310 ev 


42 km 


fluctuations of the 8th dynode current make it difficult to perform any 
more detailed analysis. Magnetometer data, however, provide a more 
quantitative comparison. 

For reasons of clarity, we shall transform the angles a and i p into 
the geomagnetic coordinate system. 

To make a comparison with the directions of the geomagnetic field 
within the radiation belt, the direction of the field was obtained from 
experimental data at three points shown in Table II. Region la occurs 
just after the recovery of the telemetry within the belt. Region lb was 
selected to have the same increment of distance as exists between III 
and la. 


Table II 


Time (U.T.) 

Region 

Radial distance (earth radii) 

17:31 

HI 

12.135 

19:42 

la 

10.877 

21:12 

lb 

9.618 


Several theoretical models indicate that in the absence of a ring 
current and close to the equatorial plane the boundary of the magneto- 
sphere is almost circular in the meridian plane on the day side of the 
earth. (Beard, I960; Beard, 1962; Beard and Jenkins, 1962; Davis and 
Beard, 1962; Mead and Beard, 1964; Spreiter and Briggs, 1962; Spreiter 
and Alksne, 1962; Spreiter and Hyett, 1963). Assuming that this holds 
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true the direction of the magnetic field was calculated just within the 
boundary of the magnetosphere at +4° geomagnetic latitude and at 10:24 
hours local time. 

The results of the field calculations and observations are given in 
Table III. 


Table IE 


Directions of the Magnetic Field 


Region 

Time 

Field(y) 

ip (deg.) 

a (deg.) 

Field Geomagnetic 
Coord. 

Lat. 

(deg.) 

Long. 

(deg). 

0 

(deg.) 

III calc. 
Ill obs-. 
Ia obs. 
Ib obs. 

17:31 

17:31 

19:41 

21:12 

61.1+0.8 

71.5+0.7 

73.5+1.1 

113.45 
1 0 5 . 0± 1 . 0 
116.3+0.7 
119.5+0.9 

126.4 

122.4+0.8 

120.4+0.6 

118.2+0.9 

86.06 

83.92 
81.11 

79.92 

139.81 
-115.50 

179.82 
166.54 

180.00 
- 75.69 
-112.53 
-107.12 


Positions of the Sun and Satellite 


Region 

Time 

Satellite Geom. Coord. 

Sun Geom. Coord. 

Lat. (deg.) 

Long, (deg.) 

Lat. (deg.) 

Long, (deg.) 

in 

17:31 

3.94 

-39.81 

10.38 

-16.08 

la 

19:41 

2.61 

-67.65 

6.82 

-48.77 

Ib 

21:12 

1.77 

-86.34 

2.79 

-71.24 


For comparison the satellite and sun positions were also included. 
A diagram of the results is given in Figure 6, where a polar projection 
is used to represent the field directions. To make possible the com- 
parison between magnetic fields as observed in different regions we 
re-define the zero magnetic meridian to lie in the plane formed by the 
earth’s magnetic polar axis and the radius vector of the satellite posi- 
\ tion. In this coordinate system the azimuthal variable of the field 

direction, <£ , is the angle whose magnitude is equal to the geomagnetic 
longitude of the field minus the geomagnetic longitude of the satellite 
\ position. 
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Figure 6-Polar plot showing the directions of the magnetic field in Regions la, lb, and III. Also 
shown is the direction of the magnetic field as calculated under assumption described in the 
text. The azimuthal variable is the difference between the geomagnetic longitudes of the field 
and of the satel life. 
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The points represent the directions of the magnetic field as ob- 
served in Regions III, la, lb, and as calculated in Region III. The circles 
about the points correspond to standard deviations due to scatter of the 
raw data points. The lack of overlap of the data is evidence for slight 
temporal variations of the magnetic field. This was confirmed by 
taking several adjacent sets of data points in both Regions la and lb. 

The results show a clustering of points with an average distance be- 
tween the points of about two degrees. 

The observed directions of the magnetic field seem to be out of 
phase with respect to the satellite by about -90° in the azimuthal di- 
rection rather than the expected 180° if the field is calculated under 
assumptions stated before. 

This difference does not seem significant. Even at 20,000 and 
30,000 kilometers field measurements disagree by about 7 degrees of 
arc with a theoretical model based on the Jensen and Cain coefficients 
and the Mead cavity field. (Jensen and Cain, 1962; Mead, 1964). This 
disagreement may be due to slight uncertainties in the calibration of 
the instruments or an error in the assumed direction of the satellite 
spin axis, or both. The problem is currently being investigated. 

The significance of these results is that in spite of several system- 
atic uncertainties in the absolute measurement of the magnetic field 
and other effects such as the change in the satellite's geomagnetic lati- 
tude from +3.94° at 17:31 U.T. to 1.77° at 21:12 U.T., and the unknown 
spatial and temporal changes in the direction of the magnetic field, the 
observed field directions in Regions III, la and lb can be fitted into a 
cone with a half angle of 3° . 

This means that the magnetic field in Region III points almost in 
the same direction as the field in the trapping region and it is, there- 
fore, reasonable to assume that Region III is just a temporary extension 
of the trapping region. 

It is also possible to put a lower limit on the speed with which the 
trapping region overtook the satellite. The initial rise in the intensity 
of the magnetic field occurred at 17:30:52.9 U.T. At 17:31:04.3 U.T. 
the detector saw 140 kev protons. During the 5.2 sec. counting period 
the detector could see protons with 90° pitch angle. They were coming 
from directions shown in Figure 7. Since these protons were trapped 
their orbit could not extend beyond the boundary of the magnetic field. 
Thus we are led to the conclusion that during the interval of 11.4 seconds 
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Figure 7-Orbits of trapped protons with 90° pitch angle which were observed in 
Region III. At the time when the trapped protons were observed the boundary of 
the magnetosphere must have been at least at a distance of 1780 km in the radial 
direction away from the satellite. The interval between the first observation of the 
increased magnetic field and the observation of 140 kev protons was 11.4 sec. 
This puts a lower limit on the velocity of advance of the magnetosphere of 156 
km/ sec. 


between the increase in the field and the appearance of 140 kev trapped 
protons the boundary must have moved through a radial distance of at 
least 17 80 km. This puts a lower limit on its radial velocity of advance 
equal to 156 km/sec. A similar analysis can be performed in trying to 
estimate the speed of recession of the magnetosphere. The time inter- 
val between the instant when ZOO kev locally mirroring protons were 
seen and the instant when the field strength dropped to about 2 5 gamma 
was 29 seconds. To preserve closed orbits of trapped protons, the 
magnetosphere must have moved in that interval at least through 2130 
km. This puts a lower limit on the radial velocity of recession equal 
to 7 5 km/ sec. 

A survey of ground magnetograms was made. No identifiable 
features could be found in the records from low latitude stations. Large 


regular pulsation events were observed in two out of approximately 
twenty high latitude magnetograms which were surveyed. At Point 
Barrow (geomagnetic latitude 68.6°N, geomagnetic longitude 241. Z) the 
field vector began rotating counter-clockwise, as seen by an observer 
looking down at the earth, at about 17:31 U.T., or 0700 local time. The 
peak to peak amplitude was about 150 gammas and the period was 4 to 
5 minutes. At Baker Take (geomagnetic latitude 73.8°N, geomagnetic 
longitude 315.2° ) the peak amplitude was again about 150 gammas and 
the period was 4 to 5 minutes, but here the wave had a roughly linear 
polarization and began about 17:27 U.T., or 11:00 local time. 

These ground records suggest that the observed rapid expansion 
and probable compression of the magnetosphere were associated with 
the excitation of a hydromagnetic disturbance in a thin shell of the outer 
magnetosphere. We do not, however, believe that the nearly simultane- 
ous observation of one event on the satellite and of an event at two out 
of twenty ground stations can be interpreted as conclusive evidence that 
the two phenomena are related. 


INTERPRETATIONS 

We have argued that at 17:31 U.T. the magnetosphere expanded past 
Explorer XII at a velocity of at least 156 km/ sec. Even if the external 
pressure dropped to zero, the plasma in the magnetosphere would not 
expand at a velocity greater than 


B 2 2c 2 
4/7mn + y - 1 


where n is the number of protons per cm 3 , c is the speed of sound, 
and we will use B equal to 60 gammas, m equal to the mass of a proton, 
and y equal to 5/3. The speed of sound in an ideal gas is given by 


c 


2 


yRT 


where R is the gas constant and j± is the molecular weight of the gas. 
Using these expressions, the expansion velocity is given by 
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+ 482T 


V 2 


1.71 x 10 6 
n 


where T is the temperature in electron volts and V is in km/sec. 

If the magnetic energy density far exceeds the thermal energy 
density, the first term will dominate and an expansion velocity of 156 
km/ sec implies a number density of less than 70 protons per cm 3 . If 
the magnetic field can be neglected, the second term dominates and 
implies a proton temperature of at least 50 ev. Several other possible 
combinations are listed in Table IV where it can be seen that either the 
number density must be less than about 100 protons per cm 3 or the 
temperature must be greater than about 15 ev near the outer edge of 
the magnetosphere. 


Table IV 


upper limit on n(cm" 3 ) 

00 

341 

174 

1 17 

88 

70 

lower limit on T(ev) 

50 

40 

30 

20 

10 

0 


After a time interval of 69 seconds a second boundary was observed. 
Several interpretations are possible. The observations could be caused 
by an expansion and subsequent compression of the entire magnetosphere, 
or a small portion of the magnetosphere; they could be caused by the 
expansion of a small portion of the magnetosphere which breaks off and 
continues to move outwards; they could be caused by a wave propagating 
in the surface of the magnetosphere. We cannot presently determine 
the cause of the observation. 


CONCLUSIONS 

On September 25, 1961 at 17:31 U.T. Explorer XII passed through 
a region of enhanced proton and electron intensity which lasted for 1 
minute and extended over an apparent distance of 83 km. At the same 
time the on board magnetometer registered an increase in the magnetic 
field. 

Analysis of the data shows that the detected protons were trapped 
in the magnetic field. At the same time their cyclotron radii were 
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much larger than the apparent width of the region. This fact points out 
the temporal character of the region. Similarities between the particle 
intensities and the strength and directions of the magnetic field in the 
region and in the outer belt suggest that the region comprises an 
impulse-like extension of the magnetosphere. 

From knowing the time of the increase and the intensity of the 
magnetic field as well as the time when the 140 kev trapped protons 
were detected, a lower limit on the outward radial velocity of the 
magnetospheric boundary can be determined to be 156 km/ sec. Simi- 
larly the lower limit on the inward radial velocity of the magnetospheric 
boundary has been determined to be 7 5 km/ sec. 

The conclusions that can be drawn from our crude calculations are 
that either the proton density near the outer edge of the magnetosphere 
is less than about 100 per cm 3 or the temperature is greater than about 
1 5 ev. 
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